function plot_SETASAnnualdiag(varname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% THIS FUNCTION PRODUCES MAPS OF THE ANNUAL MEAN AND STANDARD DEVIATION 
% OF A SELECTED SETAS MODEL VARIABLE. 
%
% INPUT     : - varname  : a SETAS variable. Available options are 'temp' 
%             or 'salt' 
% 
% EXAMPLE   :  plot_SETASAnnualdiag('temp')
%
% Author: Benedicte Pasquer, IMOS/eMII (http://imos.org.au/)
% email: benedicte.pasquer@utas.edu.au
% May 2013   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%computes setas 2yr annual mean and standard deviation
[lon,lat,setas_av,setas_std] = get_SETASAnnualMean(varname);

%FIGURE
diag2plot = {'setas_av','setas_std'};

for i = 1:length(diag2plot)
    
    plot_map(lon,lat,eval(diag2plot{i}))
    
    var4figinfo = strcat(diag2plot{i}, varname);
    
    switch var4figinfo
        case 'setas_avsalt'
            caxis([0 35]) 
            title('SETAS Sep 2010 - Sep 2012 annual mean surface salinity [ppt]')
            outputfile = 'SETAS_meanS';
        case 'setas_avtemp'
            caxis([8 16])
            title('SETAS Sep 2010 - Sep 2012 annual mean surface temperature [degC]')
            outputfile = 'SETAS_meanT';
        case 'setas_stdsalt'
             caxis([0 7])
             title('SETAS Sep 2010 - Sep 2012 salinity standard deviation [ppt]')
             outputfile = 'SETAS_stdS';
        case 'setas_stdtemp'
            caxis([0 7])
            title('SETAS Sep 2010 - Sep 2012 temperature standard deviation [degC]')
            outputfile = 'SETAS_stdT';
     end
            
    colorbar 
    saveas(gcf,outputfile,'png')
    close(gcf)
end